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ABSTRACT 


Several alternative power systems are being considered 
for small autonomous submersibles to meet long duration and 


high speed operations. The closed cycle Brayton engine 
using lithium sulfur hexafluoride as the energy source is 
one of the more promising systems. A small closed cycle 


Brayton engine has been evaluated utilizing components which 
have been tested for space applications. Modifications were 
made to optimize the unit for the undersea vehicle. 

Computer programs were developed to facilitate this process. 
It was concluded that a 2 kilowatt Brayton cycle engine 
would occupy 50 inches of length in a 21 inch diameter space 
and 36 inches in a 25 inch diameter space. 
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Chapter 1 


INTRODUCTION 


iad GENERAL REMARKS AND BACKGROUND 


Developments in the fields of navigation, controls, 
pattern recognition and artifical intelligence have made it 
possible for autonomous underwater vehicles to perform a 
wide range of tasks in the oceans. When these vehicles are 
untethered, so that power and energy have to be self- 
contained, their capabilities are seriously inhibited by the 
absence of compact power plants and energy sources. Short- 
duration and low speed missions can be accomplished with 
lead-acid or silver-Zinc batteries, but higher speed and 
longer duration missions are, at the present time, beyond 
the capabilities of existing technology. 

The general requirements of a power and energy source 
tor an effective autonomous vehicle are that it should be 
compact, safe and inexpensive. In addition, it should not 
expel exhaust products into the ocean. This latter 
requirement reduces the complexity of ballasting and 
trimming the vehicle during operation. Several studies have 
examined possible energy sources for small underwater 
vehicles [1-1] [1-2]. Recent investigations have indicated 
that a lithium sulfur hexafluoride energy system coupled to 


a small closed cycle Brayton engine provides one of the more 





et 


promising solutions to the problem. Lithium sulfur 
hexafluoride energy systems are being developed for military 
applications and small high efficiency Brayton units have 
been tested for possible space power systems for many 
thousands of hours. Fuel cells and lithium-thionyl 
batteries are serious competitors for this application [1- 
a]. 

The work described in this report continues the 
investigations of Carmichael [1-1] and Labak [1-2] on the 
lithium sulfur hexafluoride Brayton power system. The 
current work is mainly concerned with design modifications 
to the Brayton cycle system to make it suitable for 
underwater application. Much of the study was to establish 
feasible designs for the largest component of the power 
plant, the plate-fin recuperator. The recuperator can be 
configured in such a way to minimize the space occupied by 
the power plant in the vehicle. Such changes to the 
configuration have an impact on the performance of its 
Brayton cycle; and, these effects have been calculated. In 
addition, the problem of insulating the hot components has 
been studied and predictions of heat loss are provided. 
Possible system arrangements in 17 inch, 21 inch and 25 inch. 


inner working diameter vehicles are included. 
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Chapter 2 


THE POWER PLANT 


2 1 SYSTEM CYCLE 


The ability of the closed Brayton cycle power system 
to operate with good performance over a range of power 
levels and cycle conditions permits the system to be used 
for a variety of applications. The closed Brayton cycle is 
at a state of development where its performance, development 
time and cost can be accurately predicted. It can be mated 
with a wide variety of heat sources and the selection of a 
lithium sulfur hexafluoride heat source provides a very 
compact and lightweight system. 

The basic Brayton cycle consists of a compressor, an 
alternator, a turbine, a heat source and a cooler. The 
regenerative Brayton cycle, where the turbine exhaust 
temperature is quite high and exceeds the compressor outlet 
temperature, has an additional heat exchanger, referred to 
as a recuperator, that provides a higher thermal efficiency 
(fig. 2.1). A high temperature heat source, such as lithium 
sulfur hexafluoride, is required for an efficient Brayton 
cycle. 

The Mini-Brayton Rotating Unit, developed and tested 
by Garrett Airesearch [2-1], began as an outgrowth of the 


Brayton Rotating Unit for space applications. The Mini- 
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Brayton cycle is a highly developed and reliable power 
system. Considerable effort by Garrett Airesearch (for 
NASA) has been expended to demonstrate its capability. 
Garrett Airesearch’s Mini-Brayton Rotating Unit, referred to 
as the Mini-BRU, produces up to 2.4 kilowatts electrical 
Eross power and demonstrates efficiencies of 25 to 30 per 
cent. The Mini-BRU compressor and turbine are radial flow 
with a radial gap alternator mounted on a single shaft 
between them. The rotating unit operates at 52,000 rpm 
(fig. 2.2) [2-2]. The significant advantage and determining 
factor toward selection of the rotating unit is that it is 
hermetically sealed, preventing any loss of working fluid 
and providing for a clean system. The system is totally 
self-contained and bearings are lubricated by the working 
fluid, Helium-Xenon. The unit is 5.5 inches in diameter and 
13 inches in length. 

The recuperator is the largest system component in 
the Brayton cycle since it requires large amounts of heat 
transfer. A Helium-Xenon gas mixture serves as the system's 
working fluid. Because of the high heat transfer 
coefficient of this gas mixture, the size of the recuperator 
is reduced; but, it still remains the largest system 
component. Garrett Airesearch designed, built and tested a 
Mini-Brayton recuperator for NASA. Their data and 


specifications are used throughout this report [2-3]. 
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Fig. 2.2 Diagram of a Mini-Brayton Rotating Unit 
(Mini-BRU) 





Lan SYSTEM REQUIREMENTS 


System requirements investigated by this report 
included a 2 kilowatt electrical power output, mission 
durations ranging from hours to days and a vehicle speed 
equal to approximately 8 Knots. Outer diameter vehicle 


constraints examined were 21 inches, 25 inches and 34 


inches. 
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Chapter 3: 


SYSTEM COMPONENT ANALYSIS 


Sia EVALUATION OF GARRETT AIRESEARCH DATA 


Rating and sizing problems are the most common heat 
exchanger problems [3-1]. The rating problem was analyzed 
to predict the proposed heat exchanger’s performance and 
verify Garrett Airesearch’s Mini-Brayton recuperator 
specifications. Flow arrangement, core dimensions, material 
and surface geometries, flow rates and inlet temperatures 
were manually manipulated to predict outlet temperatures, 
heat transfer rate and pressure drops. Subsequently, the 
Sizing problem was used to determine exchanger dimensions by 
assuming inlet and outlet temperatures, flow rates and 
pressure drops on each side. A program was developed to 


demonstrate this objective. 


3.1.1 Rating Prokhlem Analysis 


Experimental data were analyzed to develop a 
theoretical model of the (counterflow) recuperator and to 
provide estimates of the recuperator dimensions. Detailed 
calculations are supplied in Appendix A. Surface 
geometrical properties, heat transfer coefficients and 
rates, effectiveness values, pressure drops, Reynolds number 


and fluid properties are products of a rating problem 





analysis. Since its gaseous properties are accurately 
measured, air was initially used as the working fluid. 
Recuperator dimensions and core fin configuration were used 
to derive free flow ae and total heat transfer areas. 
From these values and given mass flow rates, the hydraulic 
diameter and Reynolds number were calculated. The friction 
factor (f) and the Colburn modulus (j) values were then 
determined from the above values. Configuration properties 
and j and f curves for plate-fin surfaces, extracted from 
ref. 3-2, were compared against values and properties thus 
far calculated. For the strip-fin, plate-fin surface, all 
values compared favorably (fig. A-1). 

To verify the accuracy of the heat transfer 
equations, various core cold side mass flow rates were 
substituted into the series of equations and solved for core 
cold side pressure drops and cold side effectiveness values. 
These values were plotted on curves provided in the Garrett 
article: l)cold side air flow versus static pressure drop 
and 4)high pressure air flow versus cold side effectiveness. 
Calculated values and plotted values closely correlated. 

Satisfied with the comparison between manual 
calculations and article calculations, air properties were 
replaced with constant gas properties for a Helium-Xenon gas 
mixture. Additionally, end section pressure drops were 


included (Appendix A). Ducting and manifold pressure drops 


Ee, 





were purposely excluded since there was no accurate way to 
estimate their values. 

End section free flow areas and mid-streamline 
lengths were calculated using end section geometry and end 
section fin configuration. Each side, both the hot and the 
cold, consisted of two end sections - a hot end and a cold 
end. Pressure drops for both sides, less ducting and 
manifold pressure drops, were computed and compared to that 
predicted by Garrett Airesearch. There was good agreement 


with all values. 


3.1.2 Recuperateor Program/Siaing Problem Analysis 


Construction type and basic surface geometries had 
been selected and verified as a result of the rating problem 
analysis. Satisfied with the solution presented by the 
rating problem analysis, a program was developed to 
calculate both the hot and cold side pressure drops and 
optimally size the heat exchanger (Appendix B). Equations 
for j and f values, gas properties and material properties 
were developed. The j and f curves for the strip-fin, 
plate-fin surface were linearized for the Reynolds number 


and may be written as: 


0.008 + 6.0/Reynolds number (2321) 


Cu. 
ty 


Fh 
it 


0.03425 + 23.75/Reynolds number Core.) 
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The gas properties of helium and xenon were obtained from 
ref. 3-3. The properties for a Helium-Xenon gas mixture 
were calculated using methods from ref. 3-4 (viscosity) and 
ref. 3-5 (conductivity). The properties are temperature 
dependent and linear equations for viscosity and 


conductivity were developed as follows: 


viscosity = (@.@@0438 * T) + 1.5348 (Gaees) 
T - expressed in degrees C 
viscosity - lb/ft-sec x 19-5 
conductivity = (@.9@@038626 * T) + .916362 (3.4) 
T - degrees C 


conductivity - BTU/hr-ft-deg F 


There was some concern about the accuracy of the mixtures 
conductivity during the analysis of the Garrett Airesearch 
fomwe. since the conductivity is only used to determine the 
Prandtl number in the analysis, it was decided to 
arbitrarily increase the value of the mixture’s Prandtl 
number by a factor of 1.59. This kept all values consistent 
with Garrett Airesearch data. The recuperator and fins are 
fabricated from Hastelloy X. This material exhibits 
exceptional strength and oxidation resistance up to 2002 
degrees F. Its melting range is 2300 to 2470 degrees F. 
Its thermal conductivity was found in ref. 3-6, and may be 


written linearly as: 





Zod 


Gemduetivity = (9ae@06 * T) + 5.038 (3d 55) 
T - degress F 


conductivity - BTU/sec-ft-deg F 


The recuperator end sections are triangular shaped 
(fig. A-2) and are crossflow. To estimate their heat 
transfer values, a crossflow correction factor must be 
determined; however, this factor has only be found for flow 
that enters and exits in a parallel fashion. End section 
flow does not follow this pattern, and it is difficult to 
model. Since the heat transfer in the end sections was 
small, it was decided to neglect heat transfer in that 
region. A more conservative overall heat transfer value 
will result and essentially account for any loss caused by 
manufacturing tolerances and flow maldistribution. The 
program was designed to compute recuperator dimensions based 
on core and end section pressure drops which together 
comprise 78% of the total recuperator pressure drop. 

An initial estimate of the core cold side gas 
velocity rate, G, was determined in a loop. Core dimensions 
were determined for the side having the more strigent 
pressure drop [3-1]. Therefore, the core cold side pressure — 
drop and core dimensions immediately followed. The core hot 
side gas velocity rate was determined next from the cold 


core free flow area, assuming that the core cold and hot 
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sides possess equal free flow areas. Its pressure drop 
followed. Calculation of end section pressure drops for both 
cold and hot sides ensues. Pressure corrections to account 
for the flow acceleration, exit and entrance effects were 
determined for each pressure drop calculated using equations 
extracted from ref. 3-1. Entrance and exit effects are 
relatively small and are of opposite signs. The flow 
acceleration factor is generally less than 10% of the core 
pressure drop. All pressure drops were summed and a re- 
calculation of core cold side heat transfer values were 
determined to re-enter the main loop of the program. Each 
loop was set to a specific number of iterations discovered 
by trial and error. 

Once core dimensions were determined, additional 
length, width and height were added. These additional 
dimensions, determined from an initial trial and error run 
of the program using Garrett Airesearch pressure drops for 
calibration, accounted for end sections, manifolds and 
ducting. To parallel program and Garrett findings, only 
core cold side parameters and equations were initially 
exercised. Hot side parameters and equations were added 
when both program and article results similarly compared. 
Compact heat exchangers for closed cycle gas turbines do not 
normally have fouling problems; therefore, the fouling 
resistances were neglected. Wall resistances were found to 


be small and were also neglected. 





To ensure that the shape of the recuperator be 
maintained similarly to that of Garrett Airesearch’s, ratios 
of core width to core height and free flow area to face area 
were implemented in the program. These ratios further aided 
the determination of core dimensions. Additionally in end 
section calculations, geometry of the ends was scaled from 
the values used by Garrett Airesearch. 

Initial evaluation of the program, complete with both 
hot and cold side calculations, was performed with Garrett 
Airesearch proposed parameters and pressure drops. 

Satisfied with the program’s performance, a recuperator 


analysis was conducted. 
3a 2 PERFORMANCE PROGRAM 


A program to determine cycle efficiency, overall 
efficiency, system temperatures and system pressures was 
developed (Appendix C). This program provided input to the 
recuperator program. 

System parameters implemented were obtained from ref. 
2-1. The working fluid was a Helium-Xenon gas mixture and 
was treated as a perfect gas with constant specific heats, 
identical to the recuperator gas mixture. A two per cent 
bleed was accounted for in the computation of the 
recuperator low pressure temperature. Bearing and windage 
losses, which were density dependent, and alternator losses 


were accounted for in the determination of various system 





24 


temperatures. Detailed explanation of this program is found 


in Appendix C. 


ae, REMAINING SYSTEM COMPONENTS 


With the largest component essentially sized, 
attention was directed to the sizing of the combustor and 
the cooler and combustor heat exchangers. 

The combustor used in this system was evaluated and 
sized in ref. 1-2. It accommodates the system and vehicle 
requirements. The combustor is cylindrical in shape with a 
diameter equal to 11.2 inches and a length equal to 13.5 
inches. The combustor is the component that operates at the 
highest temperature. The amount of insulation for the 
system will be gauged from the combustor’s heat loss. Its 
heat exchanger was similarly computed and sized in ref. 1-2. 
Calculations are re-evaluated in Appendix D using the 
linearized gas property equations noted earlier. The heat 
exchanger is an inch in width and approximately 5 inches in 
length. It is oversized, in length, to accommodate any 
transient system behavior. The heat exchanger is a wavy- 
fin, plate-fin exchanger. 

The cooler heat exchanger, sea water cooling on one 
side and the Helium-Xenon gas mixture on the other, is 
designed ina similar manner to the combustor heat 
exchanger. The cooler is approximately an inch in width and 


its length varies with the vehicle diameter. Lengths range 





from approximately 16 inches to 12 inches. It is also a 
wavy-fin, plate-fin exchanger. The cooler will be hull 
mounted to save space, to alleviate the need for a pump to 
circulate the sea water and for improved maintenance. 
Cooler calculations can be found in Appendix D. 

The cooler and the combustor heat exchangers are 
fabricated from stainless steel; however, a more thermal 
conductive material, such as Hastelloy X, would improve the 
overall heat transfer coefficients and reduce component 
sizes. The reduction in these heat exchanger sizes would 


not impact the length of the vehicle powering compartment. 
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Chapter 4 


CYCLE ARRANGEMENTS 


Be) DESIGN CONSIDERATIONS 


Satisfied with the performance of the recuperator and 
performance programs, an analysis was conducted to determine 
the most compact and feasible recuperator dimensions, there- 
fore the most compact and feasible total design, while still 
maintaining overall system efficiency. The best approach to 
this task was to examine a few design considerations that 
would produce the best system arrangements. Selecting a 
specific range of pressure drop values, for use in the 
analysis, aided not only recuperator design, but system 
design as well. Additionally, to maintain system efficiency 
and to provide a viable solution to heat related problems, 
high temperature insulation must also be defined as a 


Significant design consideration in the overall analysis. 


4.1.1 Pressure Drops. and Recuperater Dimensions 


The need for the recuperator to be shorter than the 
Garrett Airesearch recuperator resulted in a range of 
recuperator pressure drops equal to or, more likely, smaller | 
than those proposed by Garrett. Evaluating the total 
system, significant pressure drops occur across the heater, 


cooler and recuperator. A @.1% pressure drop for the system 
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cooler and a @.2% pressure drop for the system heater were 
assumed. The range of total pressure drop values used were 
1.0%, @.80% and 0.50% which resulted in recuperator pressure 
drops equal to 0.70%, 0.50% and 0.30%, respectively. An 
analysis was performed and the recuperator’s length, width 
and height were computed at cold side effectiveness values 
equal to 97%, 97.5% and 98%. Table 4.1 provides a summary 
of the results. 

To arrive at final core recuperator dimensions, the 
smallest core length and corresponding face area were used. 
The height and width of the core were determined from the 
face area to provide a compact arrangement in the underwater 
vehicle. Additional length, width and height were added to 
the core to provide for recuperator casing and manifolds. 


The final overall recuperator dimensions were: 


length - 19.592 inches 
height - 11.117 inches 


width - 11.881 inches 





Delp pti * Delp pr 
1.00 0.100 
1.00 Q.10P 
1.00 0.100 
0.82 0.500 
0.80 0.500 
0.80 0.500 
0.60 0. 300 
0.60 0. 300 
0.60 0. 300 
Note: 


* Delp pt is the total system 


cent. 


Delp pr is the total recuperator pressure drop in 


TABLE 4.1 


Summary..of Recuperator Analysis 


per cent. 
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Effc is the recuperator cold side effectiveness. 


Length, 


of the recuperator in inches. 


1 


.469 
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width and height are the overall dimensions 
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4.1.2 High Temperature Insulatien 


High temperature insulation is used to reduce heat 
loss and to prevent the energy system from disturbing the 
vehicle. Several possible types of insulation products that 
would satisfy the high temperature insulation needs of the 
vehicle are discussed in the following paragraphs. 
Components requiring insulation were the combustor and the 
hot ends of the recuperator and the Mini-BRU (turbine end). 

Because the combustor is at a high operating 
temperature (about 1600 degrees F), the heat transfer is 
high and the amount of insulation required becomes great. 
Garrett Airesearch recommended using the high temperature 
insulation, graphite felt. A primary advantage of this 
material is that it can be used at temperatures well above 
3000 degrees F in an inert or vacuum environment. It 
features one of the lower thermal conductivity values 
available (in industry) and it has a low heat storage 
capacity. This insulation material is also lightweight, 
durable, easy to handle, fabricate and install. Using 2 
inches of graphite felt insulation in an inert environment 
around the combustor and the hot ends of the Mini-BRU and 
recuperator resulted in 1.2554 kilowatts of heat loss. This 
produced a 20.46% loss throughout the system. Appendix £ 
contains detailed calculations. Insulation around piping 


and hot ends of key components may not require the full 2 
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inches; however, 2 inches is provided as a conservative 
estimate. An alternative insulating material considered was 
Kaowool. It is a high temperature ceramic material with 
typical uses including gas turbines, furnances and ovens. 
Kaowool’s raw material is kaolin or alumina-silica fire 
clay. It exhibits properties and a thermal conductivity 
Similar to that of graphite felt at 1600 degrees F, 1.1 BTU- 
in/hr-ft2-degrees F and 0.98 Btu-in/hr-ft-2-degrees F, 
respectively. Other alumina-silica compositions possess 
Similar or larger thermal conductivity values at 1600 
degrees FE. 

A vacuum insulated dewar is another possible 
insulation solution. This type of insulation would require 
a high vacuum and radiation barriers. It is a clean, rigid 
and non-demountable system; however, it i8 very expensive 
and is usually fabricated to fit, and is not an off-the- 
shelf item. Another solution, offered to improve the 
insulation features of the system, was to allow a small air 
gap between hull and insulation to further reduce heat 
transfer. Care would have to be taken to ensure that 
natural convection does not become a problem in the gap. In 
conclusion, only the resistance of the insulation material 
was considered in this analysis and the results are 
therefore conservative. The contact between the combustor 
and insulation may offer some resistance, further minimizing 


any heat loss. 
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4.2 ARRANGEMENTS OF THE BRAYTON CYCLE IN SMALL 


UNDERWATER VEHICLES 


Two inches of insulation around key components was 
allowed in the following arrangements; however, diagrams 
only show insulation around the combustor. It is assumed 
that insulation around hot ends of the concerned components 
would be present. System piping is 2 inches in diameter. 
Vehicle diameters evaluated were 21 inches, 25 inches and 34 
inches with working inner diameters equal to 17 inches, Zl 
inches and 25 inches, respectively. The Mini-BRU placed 
beneath the recuperator gave the best results for piping 
arrangements and compactness. 

For a working diameter of 17 inches, the vehicle 
cannot accommodate all system components. Figure 4.1 shows 
the system arrangement with the recuperator placed 
lengthwise along the hull. Note that the cooler heat 
exchanger will not fit below the Mini-BRU and still be 
inside the 17 inch constraint, so that this arrangement is 
unacceptable. In addition, piping arrangement would be 
difficult. Placing the recuperator lengthwise across the 
hull is not possible in this diameter since recuperator 
length is greater than the working diameter. 

Figure 4.2 shows the possible system arrangement for 
the 21 inch working diameter with the recuperator placed 


lengthwise along the hull. The arrangement is tight and 
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placement of 2 inches of insulation around piping could 
prove to be difficult. Figure 4.3 is a cross-section 
looking from the turbine aft and fig. 4.4 is a cross-section 
looking from the compressor aft. 

Another alternative for the 21 inch working diameter 
is to place the recuperator lengthwise across the hull. 
Figure 4.5 and fig. 4.6 show side views of the system 
arrangement from the turbine aspect and then from the 
compressor aspect. The arrangement appears to be feasible; 
however, looking at a cross-section of the arrangement, fig. 
4.7 with the recuperator pressed up against the hull, the 
Mini-BRU protrudes out the side of the hull. 

With a working diameter of 25 inches, both system 
arrangements appear feasible. Figure 4.8 shows a side view 
of the arrangement with the recuperator placed lengthwise 
along the hull. Figures 4.9 and 4.10 show the cross- 
sectional views. Placing the recuperator lengthwise across 
the hull results in a more compact system arrangement, 36 
inches compared to 48 inches. Figures 4.11, 4.12 and 4.13 


verify the system arrangement’s feasibility. 
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Chapter 5 


CONCLUSIONS AND RECOMMENDATIONS 


5) gal SUMMARY 


The work in this report was concerned with modifying 
the design features of a Brayton cycle engine, developed for 
space applications, for use as the power source in a small 
underwater vehicle. A majority of the study focused on the 
Brayton cycle’s largest component, the recuperator, and 
established its design features from experimental results. 

A computer program was developed to provide the dimensions 
of this recuperator. An additional computer program, based 
on data for a space unit, was completed to predict the cycle 
efficiency of a small Brayton cycle engine with the 
recuperator pressure drop and its effectiveness as 
variables. The programs were used to predict engine 
performance and recuperator dimensions for a range of cycle 
parameters. A short recuperator with an appropriate height 
and width was selected for a complete Brayton cycle 
arrangement in a minimum amount of space in vehicles having 
useable diameters of 17, 21 and 25 inches. Insulation 
materials and heat losses were considered for the hot 


components of the cycle. 





Sg CONCLUSIONS 


The most compact arrangement of the Brayton cycle 
engine using the well developed Mini-BRU unit of Garrett 
Airesearch would have the characteristics shown in Table 
5.1. The thermodynamic cycle efficiency is 34.9% and the 
@veral!] efficiency is 23.1% (taking into account heat loss 
to the ocean). The overall length of the power unit is 
expected to be 5@ inches in a 25 inch diameter vehicle and 
36 inches in a 30 inch diameter vehicle. Arrangements of 
the power plants are presented in figures 4.2 - 4.4 and 


mueares 4.11 - 4.13. 


3 RECOMMENDATIONS 


The main remaining concern is the operation of the 


combustor, the storage of the fuel and the disposal of the 


products of the reaction. There is also some concern about 


insulation and heat loss. Losses could be further reduced 
by incorporating a more suitable insulating material. The 
graphite felt insulation allows about 20% loss of heat and 


could possibly be improved. 
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TABLE 5.1 


Design Characteristics of a Mini-BRU 


Working Fluid | Helium-Xenon (83.8 molecular weight) 
Net Output 2.079 kilowatts 


Pressure Ratios 


Compressor 1.491 
Turbine 1.4821 
Lost 0.994 
Efficiencies 
*Compressor 0.784 
*Turbine 0.823 
Alternator 0.99 
Cycle Q.3486 
Overall 0.281 (accounting for heat loss) 
Effectiveness 
Recuperator @.970 (high pressure) 
Recuperator @.9507 (low pressure) 


Losses (Kilowatt) 


Bearing 0.1403 
Windage Pasi ig 
Alternator OG. 1537 
Heat 1.2554 


KThese are the polytropic efficiencies. 
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APPENDIX A 
RECUPERATOR HEAT EXCHANGER CALCULATIONS 


A.1 PROPERTIES AND ASSUMPTIONS 


1. The working fluid was a Helium-Xenon gas mixture. It 

was treated as a perfect gas with constant specific heats. 
The Helium-Xenon gas mixture was the working fluid on both 
the hot and the cold side of the recuperator. 


The properties used for the Helium-Xenon gas mixture were: 


M = molecular weight is 83.8 

u ae 0.051609 x 10(-3) lbm/ft-sec 
k = 0.02889 Btu/ft-hr-deg F 

Cp = 0.05926 Btu/lb-deg R 

Pr = 0.39 

R = 0.023704 Btu/lb-deg R 


2. The strip-fin, plate-fin surface used for the 
recuperator surface had the properties shown in figure A.1 
Lose Data obtained from the NASA report was as follows: 
Core - rectangular, counterflow 

ZO-Trins/ in 


Le - fin height (plate spacing) 0.100 in 


@ - fin thickness 0.004 in 
t - plate thickness 0.008 in 
TCIN - cold temperature in 224 deg F 


THIN - hot temperature in leeZ dee F 


ol 





a2 


Mc - cold mass flow rate 09.350 lb/sec 
Mn - hot mass flow rate @.357 lb/sec 
Ee - cold side effectiveness 0.975 

PCIN - cold pressure in 196.2 psi 
PHIN - hot pressure in (iia. PSL 


49 cold plate sandwiches 


5G hot plate sandwiches 


We - core width oes as 9! 
Hght - core height 10.9 in 
Le - core length 2O...0° An 


3. Calculations were first performed with air as the 
working fluid. Pressure drops and the cold side 
effectiveness were compared with the article’s experimental 
curves for pressure drops and effectiveness. Once satisfied 
With the accuracy of the comparison and the equations, 
calculations were performed with the Helium-Xenon gas 


mixture (Table A.1). 


A.2 CALCULATIONS 


Table A.1 contains a summary of all calculation values. 


A.2.1 Eree Flow Area and Heat Transfer Calculations 


Care 


Cold side: 


Ac - free flow area 





Ac = Area rect - area fin - area plate (A.1) 
area fin = Nf * Le * Q CAS 2) 
Nf - number of fins 

Nf = # fins/in * We * # plate sandwiches CA.) 
area plate = t * We * # plate sandwiches (A.4) 

At - total heat transfer area 

Att - fin heat transfer area 

Apt - plate heat transfer area 

At = (Att + Apt) * 2 surfaces (A.5) 
Att = NE * Le * Le (A.6) 


Apt = unfinned material * Le * # plate sandwiches 


(A.7) 

unfinned material = (1 - [(# fins per inch) * Q]) 

* We (A.8) 
4rh - hydraulic diameter 
4rh = 4 * Ac * Le/At CaS), 
G - gas velocity rate 
G = mass flow rate/Ac 
Rey = G@ * 4rh/u (A.10) 


The values of the Colburn modulus (j) and the friction 


factor (f) were interpreted from fig. A.1, 


Reynolds (Rey) number. 
ROOUl =" fie * ClnIN — TCIN)] + TCIN 
TCAVG = (TCIN + TCOUT)/2 


rho - density 


entering with the 


‘Js\aale hy 
(AmeEZ) 


O3 





rho = p/(R * TCAVG) (AZ13) 
delta p - pressure drop 
delta p = £f * G(2) * 2 * Le/4rh * gravity * rho 


GAS} 


Hot side: 


Pei — tn IN ieee OUl — TOCIN) * Mc/Mn] Caveat oy) 


The hot side calculations were determined using 
equations (A.1) through (A.15) with the appropriate hot 


Side parameters used in place of the cold side parameters. 


A.2.2 Free Flow Area and Heat Transfer Calculations 
(Ends). 


1. To achieve uniform flow distribution, the average 
frictional pressure loss in the inlet end should be equal to 
the average frictional pressure loss in the outlet end. To 
obtain this distribution throughout the recuperator, it was 
decided that symmetrical end section geometries be used [2- 
3]. End sections were triangular shaped. The resultant 


flow maldistribution and performance loss were accepted. 


2. Each side of the recuperator had a cold triangular end 
and a hot triangular end (fig. A.2). 

3. The same heat transfer surface used for the recuperator 
core was used for the ends (strip-fin, plate-fin) (fig. A- 


ee 
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The end section configuration was as follows: 


Ho fas 1n 

Lte- fin height 0.100 in 
Q@e- fin thickness @.006 in 
hot end height Yeo ein 
cold end height va, 1/5 3g 


4rhe - hydraulic diameter for the end section 


4rhe = 4 X* area/perimeter 
Cold side: 
Hot end - 
Phowemn aes LOOUT ) (A.16) 
Lst - end section streamline length 


determined by means of end section geometry 


delta p = Lst K 2 * G2 * f/gravity * rho * 


4rhe (A.17) 
G = Mc/Ac(hot end) _ (A.18) 

Cold end - 
rho = p/(R * TCIN) (A.19) 


delp p = calculated as above using Ac (cold 
end) and the appropriate streamline 


length 


Hot side: 
All hot side calculations were performed in a similar 


manner to those of the cold side, both hot and cold ends, 





26 


using the appropriate lengths, mass flow rates and 


temperatures. 


To arrive at the total cold side recuperator pressure drop, 
the cold core pressure drop and the hot end and the cold end 
pressure drops of the cold side were added together. The 
hot side recuperator pressure drop was determined in the 


same Manner. 





Senitiineareti 


BEST INTERPRETATION 


aaeere 
— oa Ht 





Fin pitch = 19.82 per in = 780 per m 

Plate spacing, b = 0.205 in = §.21 x 10°2m 

Splitter symmetrically located 

Fin length flow direction = 0.125 in = 3.175 x 10°3m 

Flow passage hydraulic diameter, 47, = 0.005049 ft = 1.537 x 10°3m 

Fin metal thickness = 0.004 in, nickel = 0.102 x 10°3?m 

Splitter metal thickness = 0.006 in = 0.152 x 10°2m 

Total heat transfer area/volume between plates, 6 = 680 ft?/ft? = 2,231 m?/m3 
Fin area (including splitter)/total area = 0.841 


Fig. A.1 Properties for a strip-fin, plate-fin surface 


o¢ 
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TABLE A.1 
Summary ef Heat Transfer Calculations 


for the Recuperator 


HOT. COLD 
CORE 
TOUT (deg F) OO. 3 1260 
Nf 5800 5684 
Ac (in2) Lt PAS 21 .O7 
At (in?) SZ. oO 37511.4 
4rh (ft) 0.00537 0.00543 
M (lb/ft4-sec) 1.8838 1.862 
Rey 196.999 IRS) B) 9 tops, 
*j 0.03861 0.0386 
Kf 0.1554 0.1555 
mmo (lb/ft?) @.43597 0.684 
delta p (psi) @.1914 Ce ee, 
ENDS 
rhe (ft) O.G0641 same 
Hot 
rho Ook 9 0.4858 
delta p 0.0288 0.0181 
Cold 
rho @.72016 DPA ed 8, 
delta p 0.0652 0.0373 


*The Colburn modulus (j) and the friction factor (f) curves 


were linearized with the Reynolds number. 





APPENDIX B 
RECUPERATOR PROGRAM 


This program was developed to facilitate the determination 
of recuperator dimensions while varying pertinent input 
data. This allows the user to view several different 
recuperators simultaneously. The best ‘fit’ can then be 
determined based on pressure drops, effectiveness values and 


system requirements. 


B.1 ASSUMPTIONS 


1. The performance program (Appendix C) provided the 
temperatures, mass flow rates, pressure drop, pressures and 
effectiveness for the recuperator program's input. Input 
used from the NASA article included the fin configurations 
2-3 |. 

2. To determine a core width and core height, a ratio of 
the article’s width and height were used. This was 
Similarly done for the end section geometry. 

3. The Helium-Xenon gas mixture was considered a perfect 
fas with constant specific heats. Gas properties were 
temperature dependent and were linearized. 

4. The Colburn modulus and the friction factor were 


Similarly linearized with the Reynolds number. 
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5. As mentioned in the text, end section heat transfer was 


neglected; however, their pressure drops were calculated. 


B.2 PROGRAM 


This program read an input file and dumped data into an 


output file. A sample input file and output file follow 


the program. 





Gra? C2 


CQ) Ge Ca Ga Ge Gl C2 Ca C) Cae a Grea 1 CI) 


C2022") C2 


C202 C2 


GZ 


RECUPERATOR MAIN PROGRAM 


REAL MASSC, MASSH, LF, NTU, NTUREQ, MEANT, NCCIN, 


*NCCOUT, J, NCHOUT, NCHIN, NCC, NCH 


OPEN (UNIT = 15, FILE = ‘INPUT’ ) 
OPEN (UNIT = 5, FILE = ‘DATA’ ) 
ReaD (i>, *) KR, PCIN, THIN, “BPEL, MASSC, MASSH, LF, 


*QONE, RECUP, PCOLD, PHOT, RH, RHE, QTWO, CP 


ALL PROGRAM TEMPERATURES ARE IN DEGREES F; MASS FLOW 
RATES, IN LB PER SEC; LENGTHS, THICKNESS, WIDTHS, 
HEIGHTS, HYDRAULIC DIAMETERS IN FEET; PRESSURES AND 
PRESSURE DROPS IN PSI EXCEPT WHERE NOTED. 

R IS THE HELIUM-XENON GAS CONSTANT IN LB-FT PER LB-DEG 
R; TCIN IS THE COLD TEMPERATURE IN AND THIN IS THE HOT 
TEMPERATURE IN; EFF IS THE COLD SIDE EFFECTIVENESS ; 
MASSC AND MASSH ARE THE COLD AND HOT MASS FLOW RATES, 
RESPECTIVELY. 

LF IS THE FIN LENGTH (FOR CORE AND END); QONE IS THE 
CORE FIN THICKNESS; RECUP IS THE TOTAL RECUPERATOR 
PRESSURE DROP (CORE, ENDS, MANIFOLD AND DUCTING; A P 
ON P VALUE) PCOLD AND PHOT ARE THE COLD AND HOT 
PRESSURES, RESPECTIVELY; RH AND RHE ARE THE HYDRAULIC 
DIAMETERS DIVIDED BY FOUR CORE AND ENDS, RESPECTIVELY ; 
QTWO IS THE END FIN THICKNESS; CP IS THE SPECIFIC HEAT 
IN BTU PER LB-DEG R. 

PREQ = 0.780 * RECUP 

CORE AND ENDS COMPRISE ONLY 78% OF THE TOTAL 
RECUPERATOR PRESSURE DROP (MANIFOLD AND DUCTING 22%). 
DELP = @.20 * PREQ * PCOLD 

THE COLD SIDE PRESSURE DROP COMPRISES ROUGHLY 20% OF 
THE HOT AND COLD CORE AND END PRESSURE DROP. 

PTOTC = DELP 

NTUREQ = EFF/(1.0 - EFF) 

NTUREQ IS THE NUMBER OF HEAT TRANSFER UNITS TO EXPECT. 
HA = 2.0 * NTUREQ@ * CP * MASSC 

XJC = 6.038 

XFC = @.15 

HA IS THE HEAT TRANSFER COEFFICIENT TIMES THE TOTAL 
HEAT TRANSFER AREA AVAILABLE IN BTU PER SEC-DEG R; 

HA, XJC AND XFC ARE INITIALIZED HERE PRIOR TO THE 
LOOP. 

CALL TEMP(TCIN, THIN, THOUT, TCOUT, EFF, MEANT, MASSC, 


*MASSH, TCAVG, THAVG, EFFH) 


STARTING THE LOOP THAT REFINES THE PRESSURE DROP 


DO 106 K = 1, 12 

EP=eECOED -S8hROTC/ 2. 0 

PTOTC/2 ACCOUNTS FOR THE PRESSURE DROP FROM ENTRANCE 
TO EXIT. 
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Cymer rRer (cor, Vio, FR, KEIN, TCAVG, COND, DENS, P, R) 


THIS STARTS THE SEARCH FOR G ON THE COLD SIDE IN THE 
CORE. 


DO 180 I =1, 5 

CALL FLOW(CP, HA, G, XJC, XFC, DENS, PR, DELP, MASSC) 
CALL FRIC(RH, VI5, G, XJC, XFC) 

CALL HTRANS(EC, HC, NCC, XJC, Ghee, ek, GONE, LF, 
*KKFIN) 

CALL AREA(Z, AC, AT, DENS, XFC, G, RH, MASSC, W, AFR, 
*DELP, HGHT) 


19 CONTINUE 


CALL RECAL(XFC, Z, G, RH, DENS, PCCAL) 
CALL ETA(EHAC, EC, HC, AT) 


DETERMINING G FOR THE HOT SIDE IN THE CORE 


GH = MASSH/AC 

reer HO! om Oln/ 2 .@ 

Corie eROm@emmeyio, PR, KEIN, THAVG, COND, DENS, P, R) 
CAbby Henney lo, GH, J, EF) 

Cnet HibaNoaee, HH, NCH, J, GH, CP, PR, 

*QONE, LF, KFIN) 

CALL RECAL(F, Z, GH, RH, DENS, PHCAL) 

CALL ETA(EHAH, EH, HH, AT) 


COLD? Lwa=—--—---—= IN 


Cobh PROeG@ereeevios ER, KFIN, TCIN, COND, DENSCI, 
*XPCOLD, R) 

CAL [NRMDCWeeeAC, MASSC, RHE, DENSCI, VIS, PCIN, 
KGCIN, J, ATE) 

CALL HTRANS(ECIN, HCIN, NCCIN, J, GCIN, CP, PR, 
KQTWO, LF, KFIN) 

ECIN = (@.8633 + (@.1 * NCCIN))/@.1633 

ECIN IS THE SURFACE EFFECTIVENESS IN THE END FOR THE 
COLD SIDE; IT IS A FUNCTION OF FIN SPACING AND 
THICKNESS, IN ADDITION TO FIN EFFICIENCY. 


SOILD) (2!) == a OUT 


ees FCOLD = PIOTC 

CALL PROP(CP, VIS, PR, KFIN, TCOUT, COND, DENSCO, P, 
KR ) 

CALL OUTEND(W, AC, MASSC, RHE, PCOUT, DENSCO, VIS5, 
*GCOUT, J, ATE) 

CALL HTRANS(ECOUT, HCOUT, NCCOUT, J, GCOUT, CP, PR, 
*QTWO, LE, KFIN) 

ECOUT = (@.9633 + (@.1 * NCCOUT) )/@.1633 
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COLD PRESSURE CORRECTION 
CALL CORR(PCORRC, G, GCIN, GCOUT, DENSCO, DENSCI) 
WO Th SO SS OUT 


rae hOm — PLolTa 

Chee her Gee y Woe eh, KEIN, THOUT, COND, DENSHO, P, 
*R ) 

CALL INEND(W, AC, MASSH, RHE, DENSHO, VIS, PHOUT, 
rGHOUT, J, AIL) 

CALL HTRANS(EHOUT, HHOUT, NCHOUT, J, GHOUT, CP, PR, 
*KQTWO, LE, KEIN) 

Boul = (0.0633 (0.1 * NCHOUT))/0.1633 


OMe NW a= IN 


Cale PROD Gree vinowe rk, KEIN, THIN, COND, DENSHI, PHOT, 


*KR ) 

CALL OUTEND(W, AC, MASSH, RHE, PHIN, DENSHI, VIS, 
KGHIN, J, ATE) 

Gabe HIRANS( ERIN, HHIN, NCHIN, J, GHIN, CP, PR, 
KQTWO, LF, KFIN) 

EHIN = (0.0633 + (9.1 * NCHIN))/@0.1633 


HOT PRESSURE CORRECTION 

CALL CORR(PCORRH, GH, GHIN, GHOUT, DENSHO, DENSHI ) 
TOTAL PRESSURE SUMMATION 

Caliber hio trang, PIOTAL, PCCAL, PHCAL, PHIN, PHOUT, 
*PCIN, PCOUT, NTU, DELP, AUINV, EHAC, EHAH, 


*KEHATC, EHATH, CP, MASSC, HACAL, EC, HA, 
eco eerie een, FCORRH, NIUREQ, PIOTC, PTOTH) 
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100 CONTINUE 

an FRC CUlee nol, BEE, BRREH, PCCAL, PHCAL, Z, W, 
‘hem, FON, FOOUl, PHaaN, PHOUT, PCORRC, PCORRH, PTOTC, 
XEOUH, PTOT) 

CLOSE (15) 

CLOSE (95) 

o TOP 

END 
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RECUPERATOR SUBROUTINES 


SUMO UlMENE PROrP(Cr, VIS, PR, KFIN, [, COND, DENS, P, 


*R ) 


REAL KFIN 
DENS = (P * 144)/(R * (T + 46@)) 
VIS = ((@.@@438 * (((T + 40.8) * 5.0/9.8) - 49.@)) + 


*1.5348) * @.@0001 


COND = (©.@0@038626 * (((T + 40) * 5/9) - 40)) + 


*O.016362 


PR = 1.59 * ((VIS * CP)/COND) * 3690.2 

KFIN = ((@.0@06 * T) + 5.9038)/3600.2 

DENSITY IS EXPRESSED IN LB PER FOOTS, DENS. 
VISCOSITY IS EXPRESSED IN LB PER FT-SEC, VIS. 
GAS CONDUCTIVITY IS EXPRESSED IN BTU PER HR-FT-DEG F, 
COND. 

FIN MATERIAL CONDUCTIVITY IS EXPRESSED IN 
BTU PER SEC-FT-DEG F, KFIN. 

PR IS THE PRANDTL NUMBER. 

RETURN 

END 


SUBROUTINE FRIC(RH, VIS, G, Jd, F) 

REAL J 

REY = (4.0 * RH * G)/VIS 

J = 0.008 + 6.0/REY 

Feo 094 Zoet 2d. O/ KEY 

REY IS THE REYNOLDS NUMBER, J IS THE COLBURN MODULUS 
AND F IS THE FRICTION FACTOR. 

RETURN 

END 


SUBROUTINE TEMP(TCIN, THIN, THOUT, TCOUT, EFF, MEANT, 


*KMASSC, MASSH, TCAVG, THAVG, EFFH) 


REAL MEANT, MASSC, MASSH 

TCOUT = TCIN + EFF * (THIN - TCIN) 

THOUT = THIN - (MASSC/MASSH) * (TCOUT - TCIN) 
EFFH = (THIN -THOUT)/(THIN -TCIN) 

ani = ( (fHOUT —~ TCIN) - (THIN -TCOUT))/ 


*LOG((THOUT - TCIN)/(THIN - TCOUT) )) 


TCAVG = (TCIN + TCOW?)/2.@ 

THAVG = (THIN + THOUT)/2.@ 

ALL TEMPERATURES ARE EXPRESSED IN DEG F. 

TCOUT IS THE COLD TEMPERATURE OUT; THOUT IS THE HOT 
TEMPERATURE OUT; MEANT IS THE LOG MEAN TEMPERATURE 
DIFFERENCE; TCAVG IS THE COLD SIDE TEMPERATURE 
AVERAGE; THAVG IS THE HOT SIDE TEMPERATURE AVERAGE AND 


oye 
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EFFH IS THE HOT SIDE EFFECTIVENESS. 
RETURN 
END 


SUBROUTINE FLOW(CP, HA, G, J, F, DENS, PR, DELP, 
*MASSC) 

REAL MASSC, J 

A = (32.2 * 2.0 * DELP * 144 * DENS * J * CP * MASSC) 
B (F * HA * (PRKX (2.0/3.@0))) 

G SQRT (A/B) 

G IS EXPRESSED IN LB PER SEC-FOOT2. 

G IS THE GAS VELOCITY RATE. 

RETURN 

END 


Sono iteetinoce, H, NC, J, G, CP, PR, Q, LF, 
*KEIN ) 

REAL NC, LF, KFIN, J 

H = (J * CP * G)/(PR** (2.0/3.0) ) 

BETAC = SQRT((H * 2.0)/(KFIN * Q)) 

NC = (TANH(BETAC * LF/2.0))/(BETAC * LF/2.@) 

E = (0.05 + (@.1 * NC))/@.15 

E STANDS FOR SURFACE EFFECTIVENESS. 

NC STANDS FOR FIN EFFICIENCY. 

H IS THE LOCAL HEAT TRANSFER COEFFICIENT EXPRESSED IN 
BPUSr theta oPe-DEG R. 

BETAC IS EXPRESSED IN PER FOOT. 

RETURN 

END 


SUBROUTINE AREA(Z, AC, AT, DENS, F, G, RH, MASSC, 
*W, AFR, DELP, HGHT) 

REAL: MASSC 

Z= (2.0 * 32.2 * DENS * RH * DELP * 144)/(F * (G 
KKK2.0) ) 

AC = MASSC/G 

AT = (AC * Z)/RH 

CONST = W/H 

CONST = 0.532 

THIS CONSTANT IS TO MAINTAIN THE WIDTH AND HEIGHT 
RATIO AS IN THE NASA ARTICLE. 

SIGMA = AC/AFR 

SIGMA = 0.428 

AFR = AC/@.428 

W = SQRT(AFR * 0.532) 

HGHT = AFR/W 

AFR IS THE FRONTAL AREA; AT IS THE TOTAL HEAT TRANSFER 
AREA; W IS THE RECUPERATOR WIDTH; HGHT IS ITS HEIGHT 
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AND Z ITS LENGTH. 

AFR, AT, W, Z, HGHT ARE ALL IN FEET OR FEET2. 
RETURN 

END 


SUBROUTINE RECAL(F, Z, G, RH, DENS, PCAL) 
PCAL = (F * Z * (Gk* 2.0))/(2.0 * 32.2 * DENS * RH x 


*144) 


PCAL IS EXPRESSED IN PSI, IT IS THE PRESSURE RE- 
CALCULATED, REFINED TO ARRIVE AT A TOTAL PRESSURE. 
RETURN 

END 


THE FOLLOWING SUBROUTINE DEALS WITH ONE END OF THE 
RECUPERATOR. 


SUBROUTINE OUTEND(W, AC, MASS, RHE, P, DENS, 


*VIS, G, J, ATE) 


REAL MASS, J 

CONOUT = ZOUT/W 

CONOUT = @.3017 | 

THE ABOVE CONSTANT MAINTAINS THE RATIO OF THE END‘S 
LENGTH TO THE CORE’S WIDTH AS IN THE NASA ARTICLE. 
ZOUT = W * 0.3017 

CONTWO = ACEND/ACCORE 

CONTWO = @.6143 

THE ABOVE CONSTANT MAINTAINS THE RATIO OF THE END‘S 
FREE FLOW AREA TO THE CORE’S FREE FLOW AREA AS IN THE 
NASA ARTICLE. 

ACEND = @.6143 * AC 

ATE = (ACEND * ZOUT)/RHE 

MASS/(AC * @.6143) 

= (4.0 * RHE * G)/VIS 

@.03425 + 23.75/REY 

@.008 + 6.0/REY 

(ZOUT * F * (Gk* 2.0))/(32.2 * DENS * RHE * 144 x 


0 


ro Gi ty G2 
mom om ew 


*K2.0) 


ZOUT IS THE END LENGTH; ACEND IS THE END FREE FLOW 
AREA; ATE IS THE END TOTAL HEAT TRANSFER AREA AND G, 
REY, J, F, AND P ARE DEFINED AS BEFORE BUT APPLY 
TO THE “END* HERE. 

UNITS ARE DEFINED AS ABOVE IN PREVIOUS ROUTINES. 
RETURN 

END 
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SUBROUTINE ETA (EHA, E, H, AT) 

EHA = E * H * AT 

EHA IS EXPRESSED IN BTU PER SEC-DEG R. 

IT IS THE SURFACE EFFECTIVENESS TIMES THE HEAT 
TRANSFER COEFFICIENT TIMES THE TOTAL HEAT TRANSFER 
AREA. 

RETURN 

END 


THE FOLLOWING SUBROUTINE DEALS WITH THE OTHER END OF 
THE RECUPERATOR. 


SUBROUTINE INEND(W, AC, MASS, RHE, DENS, VIS, 
KE, G, J; ATE) 

REAL MASS, J 

CONIN = ZIN/W 

CONIN = 0.2543 

THE CONSTANT ABOVE IS THE RATIO OF THE END LENGTH TO 
THE CORE WIDTH AS IN THE NASA ARTICLE. 

ZIN = @.2543 * W 

CONTHR = ACEND/ACCORE 

CONTHR = 9.1829 

THE ABOVE CONSTANT IS THE RATIO OF THE END’S FREE 
FLOW AREA TO THE CORE’S FREE FLOW AREA AS IN THE NASA 
ARTICLE. 

ACEND = @.1829 * AC 

ATE = (ACEND * ZIN)/RHE 

G = MASS/(AC * 0.1829) 

REY = (4.0 * RHE * G)/VIS 


F = 0.903425 + 23.75/REY 

J = 0.008 + 6.0/REY 

P = (ZIN * F * (Gk* 2.0))/(32.2 * DENS * RHE *144 x 
*Z.0) 


THE PARAMETERS AND UNITS ARE THE SAME AS IN THE 
SUBROUTINE OUTEND BUT APPLY TO THE OTHER END OF THE 
RECUPERATOR. 

RETURN 

END 


THE SUBROUTINE PRES IS A SUMMATION OF THE PRESSURES; 
IT REFINES THE COLD SIDE PRESSURE DROP AND ‘HA’ VALUES 
TO ENTER THE LOOP AGAIN. 


SUBROUTINE PRES(PREQ, PTOTAL, PCCAL, PHCAL, PHIN, 
*PHOUT, PCIN, PCOUT, NTU, DELP, AUINV, EHAC, EHAH, 
*KEHATC, EHATH, CP, MASSC, HACAL, 

‘ie te, ECOGD, PHOT, PCORRC, “PCORRH, NTUREQ, PTOTC, 
*PTOTH ) 

REAL NTUREQ@, NTU, MASSC 

PTOTC = PCCAL + PCIN + PCOUT + PCORRC 
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Miia seneaib + FHIN + PHOUI + PCORRH 
EleraL = PTOTC/PCOLD + PTOTH/PHOT 
DELP = PCCAL * (PREQ/PTOTAL) 


KHATC = EHAC 
EHATH = EHAH 
AUINVY = (1/EHATC) + (1/EHATH) 


NTU = 1/(CP * MASSC * AUINV) 

HACAL = EHAC/EC 

HA = HACAL * (NTUREQ/NTU) 

PTOTC IS THE TOTAL COLD SIDE PRESSURE DROP (CORE AND 
ENDS). 

PTOTH IS THE TOTAL HOT SIDE PRESSURE DROP (CORE AND 
ENDS). 

PTOTAL IS THE TOTAL RECUPERATOR PRESSURE DROP (HOT AND 
COLD, CORE AND ENDS ONLY). 

EHATC AND EHATH ARE VALUES ONLY FOR THE CORE (HOT AND 
COLD). 

HEAT TRANSFER FOR THE ENDS WAS NEGATED, AS STATED IN 
THE ASSUMPTIONS. 

AUINVY IS THE VALUE OF 1/AU IN SEC-DEG R PER BTU. 
HACAL IS A REFINING OF THE HA VALUE AS IS DELP. 
RETURN 

END 


SUBROUTINE CORR(PCORR, GAVG, GIN, GOUT, DENSOT, 


*XDENSIN) 


REAL KE, KC 

KC = @.33 

KE = 8.31 

SIGEND = @.433 

SIGEND IS SIGMA OF THE ENDS 

FREE FLOW AREA = W/G 

SIGMA = FREE FLOW AREA/FRONTAL AREA 

PFLO = ((GAVG** 2.0)/32.2) * ((1/DENSOT) - (1/DENSIN)) 
PENT = ((GIN** 2.0)/(2 * 32.2)) * (1 - (SIGEND** 2.0) 


*+ KC)/DENSIN 


PEXT = ((GOUTK* 2.0)/(2 * 32.2)) * (1 - (SIGEND** 2.0) 


*- KE) /DENSOT 


PCORR = (PRLO + PHNT - PEXT)/144 

PFLO, PENT AND PEXT ARE PRESSURE CORRECTIONS FOR THE 
FLOW, THE ENTRANCE AND THE EXIT, RESPECTIVELY. 

PCORR IS A SUMMATION OF THESE CORRECTIONS IN PSI. 
RETURN 

END 


THE PRINT SUBROUTINE NOT ONLY PRINTS OUT THE OUTOUT 
FILE BUT ALSO ENSURES THAT PERTINENT VALUES ARE IN THE 
CORRECT UNITS. 





SUBROUTINE PRINT(PCOBD, PHOT, EFF, EFFH, PCCAL, PHCAL, 
*Z, W, HGHT, PCIN, PCOUT, PHIN, PHOUT, PCORRC, PCORRH, 
heuer PLOTH, PTOT) 

PTOTC = (PCCAL + PCIN + PCOUT + PCORRC) * 10@/PCOLD 

PTOTH = (PHCAL + PHIN + PHOUT + PCORRH) * 100/PHOT 

EROLe= (PLOre + PLORy 100 

POeCAL PCCAL * 100/PCOLD 

PHCAL PHCAL * 100/PHOT 

rem = POINe* 16¢7ePCOLD 

PCOUT = PCOUT * 100/PCOLD 

PHIN = PHIN X* 100/PHOT 

PHOUT = PHOUT * 10@/PHOT 

PCORRH = PCORRH * 10@0/PHOT 

PCORRG = PCORRC * 1G07PCOLD 

jae L - <a? 

HGHT = HGHT x 12 

W-= W X* 12 

OZ = 3,-+ baz 

OHGHT = HGHT + 1.1 

OW = W + 2.9 

OZ, OW AND OHGHT ARE THE OVERALL DIMENSIONS WHEREAS Z, 

W AND HGHT ARE THE CORE DIMENSIONS, ALL IN FEET. 

Wile (5, 165) 

65 FORMAT (/5X, “HEAT EXCHANGER DATA’ //5X, “TOTAL PRESSURE 
*XDROP’ 4X, “EFFECT. COLD HOt) 
Wigktneo, fo) Flol, EFF, ERFH 
TEMMPORMAT (Gk, SRIG.5, 12X, 2F12.4) 
WRITE (5, 738) 
78 FORMAT (/5X “PRESSURE DROP, PER CENT’ 
*K/8X,° COLD CORE HOT CORRECT. 
KTOTAL’ ) 
Nieshh wGoeeeo? )ereunN. ECCAL, PCOUT, PCORRC, PTOTC 
WRITE (5, 97) PHIN, PHCAL, PHOUT, PCORRH, PTOTH 
OF PORMAT (5X, 5h11.5) 

WRITE @5 5105) 

WRITE (5, 106) Z, W, HGHT 

WRITE (5, 106) OZ, OW, OHGHT 

105 FORMAT (/11X, “DIMENSIONS INCHES CORE AND HEAT EX- 


*KCHANGER/ 16X “LENGTH WIDTH HEIGHT * ) 
106 FORMAT (16X, 3G10.5) 
RETURN 


END 
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Sample Input File for the Recuperator Program 


ema 2iOl4, 1335.6, 2.975, 02.3328, 8.3396, 0.0083, 
@.00033, 0.00696, 99.871, 67.4, 0.081313, O@.00154, BO.0GO50, 


0.05946 


where 
R - 18.442 1ib-ft/lb-deg R 
XTCIN - 218.4 deg F 
XKTHIN - 1335.6 deg F 
-KEFF - @.975 
KMASSC - 0.3328 1b/sec 
KXMASSH - 0.3396 lb/sec 
LF - @©.0083 feet 
QONE - 0.00033 feet 
*KRECUP - @.00696 
*PCOLD - 99.871 psi 
*KPHOT - 67.4 psi 
RH - @.9@01313 feet 
RHE - 0.00154 feet 
QTWO - 0.00050 feet 
CP - @.@5946 BTU/1b-deg R 


KInput provided by the performance program 


al 
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Sample Output File for the Recuperator Program 


HEAT EXCHANGER DATA 


TOTAL PRESSURE DROP KFFECT. COLD HOT 
.O9543 .9750 Fe) OF oye 
PRESSURE DROP, PER CENT 
COLD CORE HOT CORRECT TOP RG 
.03363 10555 .O2019 01210 .17147 
.G4664 . 29060 8341 SEO OIZ3 -ST141 


DIMENSIONS INCHES CORE AND HEAT EXCHANGER 
LENGTH WIDTH HEIGHT 
23.242 SOLO 19.391 
29.442 8.4282 11.491 





APPENDIX C 


PERFORMANCE PROGRAM 


In addition to providing the necessary recuperator program 
input, the performance program determined the overall 


efficiency and cycle efficiency. 
C.1 ASSUMPTIONS 


1. Bearing and windage losses were considered to be density 


dependent. 


2. Design data used was from a NASA article on the design 


and fabrication of the Mini-BRU [2-1]. 
C.2 PROGRAM 


This program read an input file and dumped data into an 
output file. A sample input file and output file follow the 


progran. 


fe 
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PERFORMANCE MAIN PROGRAM 


REAL MRAT, MCOMP, MTURB 

OPEN (UNIT = 10, FILE = IN’) 

OPEN (UNIT = 15, FILE = ‘OUT’) 

Rebelo. +) Tis fo, Fi, GFR, MRAT, PNET, EFFC, DELP, 


*XRECUP 


T1 IS THE COMPRESSOR INLET TEMPERATURE. 
T6 IS THE TURBINE INLET TEMPERATURE. 
P1 IS THE COMPRESSOR INLET PRESSURE. 
CPR IS THE COMPRESSOR PRESSURE RATIO. 


MRAT IS THE RATIO OF THE TURBINE MASS FLOW RATE TO THE 


COMPRESSOR MASS FLOW RATE. 

PNET IS THE NET POWER OUTPUT. 

EFFC IS THE RECUPERATOR COLD SIDE EFFECTIVENESS. 
DELP IS THE TOTAL SYSTEM PRESSURE DROP. 

RECUP IS THE TOTAL RECUPERATOR PRESSURE DROP. 
TEMPERATURES ARE IN DEGREES R, PRESSURES ARE IN PSI, 
POWER IS IN KW, DELP IS IN PERCENT, MASS FLOW RATES 
ARE IN LB PER SECOND, RECUP IS IN PERCENT. 

TPR IS THE TURBINE PRESSURE RATIO. 

TPR = CPR * (1.0 - (DELP/1900.@)) 

PGROSS = PNET/@.9724 

@.9724 IS THE NET CONTROL ANALYSIS. 

P2 = Pl * CPR 

MCOMP = 90.4 

P2 IS THE COMPRESSOR OUTLET PRESSURE IN PSI. 

MCOMP EQUALS @.4 TO SET AN INITIAL BALLPARK VALUE TO 
ENTER THE LOOP, IT IS THE COMPRESSOR MASS FLOW RATE. 


STARTING LOOP TO REFINE MCOMP 


DO 100 I= 1, 12 

MTURB = MCOMP * MRAT 

MTURB IS THE TURBINE MASS FLOW RATE. 

Cali COMES (TIe@ER,sBrC, TZ) 

CAO LODE 1 Omen Wil, 7 ) 

Comte oGoom nz yeice li, GROSS, MCOMP, MTURB, T3, T4, 


*T8 ) 


T9 = MRAT * T8 + (1.0 - MRAT) * T4 
T9 IS THE LOW PRESSURE INLET TEMPERATURE TO THE 
RECUPERATOR, IT IS INFLUENCED BY A 2% BLEED (THIN). 


CALL REGEN (T4, T9, EFFC, MRAT, T5, T1i®@, EFFH) 
CADE OhwGliera TS, T6, TS, PGROSS, PNET, MRAT, 


*MCOMP, EFFCYC, EFFIC, P6, P2, P9, CPR, RECUP) 


THE PERFOR ROUTINE REFINES THE VALUES OF P2 AND MCOMP 
TO BE CARRIED OVER AT THE START OF THE LOOP. 


190 CONTINUE 
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eee i leizeto. 24, 75, 16, T7, T8, TY, T1e, 
Pie teetonCyY eC snp hto, DELP, MIURB, EFFC, P6, P2, P9) 
CLOSE (10) 

CLOSE (15) 

MCOMP AND MTURB ARE THE HOT AND COLD MASS FLOW RATES, 
RESPECTIVELY, FOR THE RECUPERATOR. 

END 
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PERFORMANCE SUBROUTINES 


SUBROUTINE COMP (il, CPR, EPC, T2) 

PR = CPR 

EPC = 0.78397 

EPC IS THE PRE-DETERMINED POLYTROPIC EFFICIENCY FOR 
THE COMPRESSOR. 

EP = 1/EPC 

CALL PROP (PR, EP, TRAT) 

T2 = Tl * TRAT 

T2 IS THE COMPRESSOR OUTLET TEMPERATURE. 

PR STANDS FOR PRESSURE RATIO, EP STANDS FOR POLYTROPIC 
EFFICIENCY AND TRAT IS THE TEMPERATURE RATIO. 
RETURN 

END 


SUEROUMINEe URED (16, TPR; EPT, TT) 


Ea Seales 
EPT = 0.82174 
BP = EPT 


Gabin PROBS( ER, EP, TRAT) 

T7 = T6/TRAT 

BPT Is THE PRE-DETERMINED POLYTROPIC EFFICIENCY FOR 
THE TURBINE AND T7 IS THE TURBINE OUTLET TEMPERATURE. 
RETURN 

END 


SUBROUTINE REGEN (T4, T9, EFFC, MRAT, TS, T10, EFFH) 
REAL MRAT, MCOMP, MTURB 

ioe [4 + (BERG * (T9 -T4)) 

T1090 = T9 - (TS -T4) * MRAT 

Been = (95> ie Cig — T4) 

TS IS THE HIGH PRESSURE OUTLET TEMPERATURE FROM THE 
RECUPERATOR, T1@ IS THE LOW PRESSURE OUTLET TEMPER- 
ATURE FROM THE RECUPERATOR AND EFFH IS THE RECUPERATOR 
HOT SIDE EFFECTIVENESS. 

RETURN 

END 





CGC 


C20 


Pe OP Bee) 


Cire C2 


SUBROUTINE PROP (PR, EP, TRAT) 

R = 18.44 

CP = 0.05946 

CP IS IN UNITS OF BTU PER LB-DEG R, SPECIFIC HEAT. 
R IS IN UNITS OF FT-LB PER LB-DEG R, HELIUM-XENON GAS 
CONSTANT. 

GAMMA = 5.0/3.2 

EXPO = (GAMMA - 1) * EP/GAMMA 

EXPO STANDS FOR EXPONENT. 

TRAT = PR** EXPO 

RETURN 

END 


SeKOUlmwe Beem (i 2qeta, 17, PGROSS, "COMP, MTURB, T3, 
*T4, T8) 

REAL MRAT, MCOMP, MTURB 

BEARING LOSS IS @.154 KW OR @.162 BTU PER SEC. 

WINDAGE LOSS IS 0.085 KW OR @.090 BTU PER SEC. 
ALTERNATOR LOSS IS 8 PERCENT OF THE GROSS POWER. 
LOSSES ARE BASED ON DENSITY (DENSITY AT 2/DENSITY AT 


BLOSS = 0.162 * PZ * 665.0/(107.0 * T2) 
= 9.9090 * PZ * 665.0/(107.0 * T2) 

ALTL = PGROSS * 9@.98/1.955 

ONE BTU EQUALS 0.9481 KW-SEC. 


CP = 0.05946 
T3 = T2 + (BLOSS + WINDL)/(2.0 * MCOMP * CP) 
T4 = T3 + (ALTL/(MCOMP * CP)) 


Totti Cet Omoet+ KIND) /(2.0 * MTURB * CP) 

T4 IS THE HIGH PRESSURE RECUPERATOR INLET TEMPERATURE 
(TCIN). 

RETURN 

END 


UP OULUNDwEnREORer dd, td, T5, T6;, T8, PGROSS, PNET, 
Tol NCOP wen MGN Cw sEnBHIC, P6, PZ, PS, CPR, RECUP) 
REAL MRAT, MCOMP, MTURB 

CP = 9.05946 

MCOMP = PGROSS/(CP * 1.055 * ((T6 - T8) * MRAT + (T1 - 
*T4))) 

P6 = (MCOMP * (SQRT(T6)))/@.15499 

THE CONST @.15499 IS IN UNITS OF DEG R TO THE SQRT 
LBM-IN2 PER SEC-LBF, THIS IS THE VALUE OF THE MASS 
FLOW FUNCTION THAT IS BEING IMPOSED. 

P2 = P6/(@.998 - @.31 * (RECUP/1@®@) ) 

P9 = (P2 * (1/CPR)) * (1.001 + @.69 * (RECUP/1909) ) 
EFFCYC = ((T6 - T8) * MRAT + (T1 - T4))/(MRAT * (T6 - 
*T5) ) 





Gare? C2 C) GIG) lA @ 


C2 2 


fs: 


PEC o=——iitt/ (tina + MCOHP * CP * (T6 -T5) x 1.855) 
P2 AND P9 EQUATIONS DETERMINED THROUGH SYSTEM PRESSURE 
DROPS (HEATER - @.2%, COOLER - @.1%) AND THE 
ASSUMPTION: PRESSURE DROP ON THE HOT SIDE OF THE 
RECUPERATOR EQUALS 0.69 TIMES THE TOTAL RECUPERATOR 
PRESSURE DROP. . 

EFFCYC 1S THE CYCLE EFFICIENCY AND EFFIC IS THE 
OVERALL SYSTEM EFFICIENCY. 

P2 IS THE RECUPERATOR COLD PRESSURE INLET AND P9 IS 
THE RECUPERATOR HOT PRESSURE INLET (PCOLD AND PHOT). 
RETURN 

END 


SUpnewmINE PRINT (Til, Iz, T3, T4, 15, T6, T7, T8, TY, 
Tit eNOOnE eh EOCYO, BERIC, DELF, MTURB, EFFC, P6, P2, 
*PY) 

REAL MRAT, MCOMP, MTURB 

WRITE (15, 150) 

WRITE (15, 200) EFFCYC, MCOMP 

WRITE (15, 200) EFFIC, MTURB 

WeeebE (15, 225) 

Wovemoto, 42400) Tl, T2, T3, T4, TS 

wet owecoe) 16, If, 18, T9, T12d 

Hrelerec lo, 26d) 

VivVQttwieroye 265) KMBG, DELP, P6, P2, PY 

150 FORMAT(10X, “GAS TURBINE DATA’//10X “EFFICI (6 SPACES) 
*KMASS° 1X° FLOW’ ) 
200 FORMAT(2X, 2F13.4) 
moo BORMAT(C/AGx, TEMPERATURES DEG R T1 TO T1d) 
265 FORMAT(1@0X, 2F8.4, 4X, F8.4, 5X, F8.4, 5X, F8.4) 
250 FORMAT(10X, 5F8.1) 
260 FORMAT(/10X, ~ EFFEC DELP Ins In {@alal 
*KSPACES) P9° ) 
RETURN 
END 





Sample Input File for the Performance Program 


eee OOo eels?! , ©.96, 2.079, 0.975, 0.996, 


@.696 

where 

T1 - 542 deg R 
T6 - 2060 deg R 
mie — (leaf psi 
Cee - 1.491 
MRAT - 0.98 
pee - 2.079 KW 
ba - 0.975 
DELP - 2.996% 
RECUP - 9.696% 





Sample Output File for the Performance Program 


GAS TURBINE DATA 


HREICI MASS FLOW 
. 35983 . 3396 
. 3407 . 3328 


TEMPERATURES DEG R Ti TO T10 
542.0 664.5 670.4 678.4 1767.7 
2060.0 1812.5 1818.4 1795.6 it Zier 


EFEC DELP P6 eZ see) 


. 9750 . 9962 99.4554 99.8706 67.3710 
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APPENDIX D 


COMBUSTOR AND COOLER HEAT EXCHANGER CALCULATIONS 


D.1 PROPERTIES AND ASSUMPTIONS (COMBUSTOR HEAT 


EXCHANGER) 


1. Reference 1-2, Appendices C and D, was used as the 
basis for these calculations. The assumptions made 

there apply similarly to these calculations unless 
otherwise noted. 

2. The heat transfer required by the Brayton cycle was 
6480 W or 22,109 BTU/hr (q), the mass flow rate was 
@.3258 lb/sec and the log mean temperature difference was 
313.69 degrees F [1-2]. 

3. At a pressure of 97.6 psia (combustor inlet) anda 
temperature of 1299.9 degrees F (combustor inlet), the 


Helium-Xenon gas mixture properties were as follows: 


rho - density 0.4330 1b/ft3 
u - viscosity 4.620 x 10-9 l1b/ft-sec 
k - gas conductivity Q@.04357 BTU/hr-ft-deg F 
Cp - specific heat @.05946 BTU/lb-deg F 
Pr - prandtl number @. 3609 

4. The heat exchanger was a counterflow configuration 


With a wavy-fin, fin-plate surface exhibiting the 
properties shown in figure D.1 [3-2]. 


5. The combustor’s diameter was 10 inches. 


ll 





D.2 COMBUSTOR HEAT EXCHANGER CALCULATIONS 


Table D.1 contains a summary of all calculation values. 
Ac = (area between the two shells) - (cross-sectional 


area of all the fins) 


(mw * (Ro2 - Ri2) - [ENF * Fin length * fin 
thickness }) (D.1) 
nie = ) inches 


Ro = 5 + @.413 = 5.413 inches 


Nf = 2* wrk Ri * # fins per inch DZ ) 
G = M/Ac Gb 3) 
M- mass flow rate 

Rey = G@ * 4rh/u (D.4) 


The friction factor (f) and Colburn modulus (j) were 
determined from the curves in figure D.1, entering with 
the Reynolds number. 
h = j * Cp * G/Pr2/s3 CDr 5) 
h - heat transfer coefficient 
n = tanh (B * Le)/(B * Le) (D.6) 
where B = [(h * 2)/Kss * Q)]!/2 Cie 7 ) 
Le - fin length 
n - fin efficiency 
Q@ - fin thickness 


Kss - fin material conductivity 
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E = (0.05618 + (2 * @.413) * n)/(@.@5618 + 2 * G@.413) 
OID te) 
E - surface effectiveness, dependent on fin geometry 
LMTD/q = 1/(Arftt * E x h) + x/(Kss * Anft) (D.9) 
where 
x - shell thickness, 0.2 inches 
Att = L * NE * 2(sides) * Ler oD 10) 
Anftt = L *®* (w * Di - NE * Q) Cari) 
LMTD - log mean temperature difference 
q - heat transfer required by Brayton Cycle 
All values necessary to solve equation (D.9) for L 
(heat exchanger length) were calculated in equations 


meet) through (D.8). 
D.3 PROPERTIES AND ASSUMPTIONS (COOLER HEAT EXCHANGER) 


1. The cooler was a parallel flow configuration with 

the exact same heat transfer surface as the combustor 

heat exchanger. 

2. The log mean temperature difference used was 77.71 
degrees F; the mass flow rate was 0.3324 lb/sec; the shell 
thickness equaled @.188 inches (hull thickness). 

3. At 65.7 psia (cooler inlet) and 273.6 degrees F (cooler 


inlet), the Helium-Xenon gas mixture properties were: 


rho @.@6993 lb/ft 


u Zee xe EO-> 1lb/it-sec 


83 





84 


k @.@2155 BTU/hr-ft-deg F 
Cp @.@5946 BTU/lb-deg F 
Pr Q.3352 


D.4 COOLER HEAT EXCHANGER CALCULATIONS 


Table D.1 contains a summary of all calculation values. 
Cooler calculations were determined using equations 

(D.1) through (D.11). Only equation (D.9) was somewhat 
modified by the addition of another term to account for 
additional resistance created by the convection of heat 


through a boundary layer in the sea water [1-2]: 


+ 1/(E * h * A) 
where 
A=2%* wt * Rv * L Chet Z) 
A - surface area of vehicle 
Rv - radius of vehicle 
KE= 1 


h 


106 BTU/hr-ft2-deg F 
h - ‘conservative’ convective coefficient 
estimate of one knot and 2 kilowatts 


power output [1.2] 


The cooler had to be designed for more than one size 


vehicle; therefore, the radii were as follows: 


1) Rv = 18.5" Ro = 9.5" Ri = 9.98" 
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Z) NG Siem Ro = 14" Rae] 1a 587: 


Seeky = 12.5" Rowe lh. >” Ri = 11.987" 


NOTE: The gas mixture properties were temperature 
dependent and were calculated using linearized 


equations. 





2° 
Pe ie 
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BEST INTERPRE TATION 
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Fin pitch = 17.8 per in = 701 perm 

Plate spacing, b = 0.413 in = 10.49 x 10°3m 

Flow passage hydraulic diameter, 47, = 0.00696 ft = 2.123 x 10°3m 

Fin metal thickness = 0.006 in, aluminum = 0.152 x 10°3m 

Total heat transfer area/volume between plates, B = 514 ft?/ft? = 1.686 m27/m3 
Fin area/tota! area = 0.892 


Note: Hydraulic diameter based on free-flow area normal to mean 
flow direction. 


Fig. D.1 Properties for a wavy-fin, plate-fin surface 
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mABLE D.1i 
Summary_of Heat. Exchanger Calculations 


for the Combustor and Cooler Heat Exchanger 


COMBUSTOR COOLER. (1). 

Ac (fit2) Zee 2 1 Zoe 9S 
Nf 299 1016 

G@ (lb/in2-sec) 0.02687 0.01508 
Rey O02 .90 712.01 

J OE One 2 0.0149 

3m 0.076 0.069 

h (BTU/hr-ft2-deg F) 24.836 14.353 

B (per ft) 87.418 818), AN) o 

n 0.331 0.4283 

E 0.374 0.4647 
Att (L * in) 461.734 6) 6s) VArA 
Anft (L * in) Zor Co oa 

Av (l * in) = Oyo ST 
LMTD/q (hr-deg F/BTU) @.9014188 0.003515 
L (in) TA OR il Ze sy 


*4.25 (oversized) 


*The combustor heat exchanger will be oversized 


approximately 5@% to account for any transient behavior 


fol= 2)... 





Ac (ft2) 
Nf 


feeenb/in2-sec) 


Rey 
J 
f 


m (BTU/nr-ft2-dege F) 


B (per ft) 
n 


K 


Att (L * in) 

Anft (L * in) 

Av (L * in) 

LMTD/q (hr-deg F/BTU) 


L (inches) 


TA Binge. 1. (cont ) 


COOLER (2) 


DIZOL 
1522 
@.@1@175 
480.41 
0.0175 
0.08 
11.374 
Se) dhe Ks, 
0.4747 
0.5982 
IASYS) = 874 
FS Ae 
a= 25 
0.003515 


10.63 


COOLER. (3). 


Ze 23 
1240 
0.@12421 
586.46 
O.O17 
0.078 
13.488 
64.422 
0.4404 
0.4760 
1924.24 
SL (AVE 
78.54 
same 


Lace 
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APPENDIX E 


INSULATION CALCULATIONS 


E.1 PROPERTIES AND ASSUMPTIONS 


1. The following table of thermal conductivity values 


(k) for graphite felt were used: 


TABLE E.1 


Thermal. Conductivities. for Graphite Felt 


k (BTU-in/hr-ft-deg F) T (deg F) 
0.8 10020 
1.1 ZOOO 
Ines 3OBO 
3.0 49000 
aya, VA DOOD 


These values were obtained from a commercial supplier 
of graphite felt. 

2. Components considered for insulation were the 
combustor and the hot ends of the recuperator and the 
Mini-BRU (turbine end). 

3. The combustor was assumed to be at a constant 
temperature of 1600 deg F. The log mean temperatures 
were calculated for the hot ends of the recuperator and 


vem eras 0) = 





E.2 CALCULATIONS 


Combustor - cylindrical in shape 
diameter - 11.2 inches length - 13.5 inches 
Area = (2 * wm *¥ R *¥ H) + 2 *® (w * R2) = 4.667 £ft2 
eH} 
R - radius H - height 
Interpolating from Table E.1, the value of "k" for 1602 
degrees F was @.980 BTU-in/hr-ft2-deg F. 
Heat loss for 1600 degrees F, 4.667 ft2 and 1 inch 


insulation equaled 


(0.986 * 1600 * 4.667 * .2931)/1 


2.145 kilowatts 
0.2931 was used for the conversion from BTU/hr to W. 
Heat loss with 2 inches and 3 inches of insulation was 


1.0724 and @.715 kilowatts, respectively. 


Mini-BRU - treated as a rectangular shape 
width - 5.5 inches height - 5.5 inches 
Area = Width * Height = @.2101 ft2 
Log mean temperature = (Ts - T7)/in((Ts - 85)/(T7 - 
85) ) 
(eZ ) 
Te was the turbine inlet temperature and T7 was the 
turbine outlet temperature. They were 1600 and 1350.1 
degrees F, respectively. The 85 degrees F was the 


assumed seawater temperature. The log mean temperature 


IG 





% AQT * ) » a ™ + . 4 = = ~s ‘ = 2 
equaled 1d8e.93 degrees § iInterpoleting rtrowe Tsble £.1 
. =~ We QV ve = ™~ Aes = . a % — P+ F . a b a) vy 
the value of “k" was @.9153 BTU-invhr-fte-des FF. Hest 


> _ ra wT os a y i ol = oC a ra Pe . re ‘ 
less for 1096.3 degrees F. OO. 2101 Ft snd 1 inch 


@.03991 and O.0261 Kilowatt, respectively. 


Recuperator - treated as a rectansular shape 


width - 11.881 inches height - 11.117 inches 
Area = @.9172 rt 


The los mean temperature was calculated using equation 


Y mall 


(KR.2) with THIN replacing Te and TCOUT replacing T?7. 
Their values were 1988.9 and 1299.9 degrees F. 


respectively. The log mean temperature was i2301.4 


=~] 
{3 
4, 
Ps 

7) 


degrees F; so for that log mean temperature, 0.91 


QO. 2878 


+4 


and 1 inch insulation, the heat loss wa 


Ki 


kilowatt. For 2 ahd 3 inches of insulation, the heat 


less was 0.1439 and 0.0950 kilowatt, respectively. 


The amount of insulation used was based on the 
degree of heat loss of the combuster (the hottest 
camponent and therefore the one with the most heat 
loss) with the varying thicknesses of insulation and 


the feasibility of system arrangement with this ameunt 


fa 





of insulation. 
So for 2 inches of insulation, the total heat loss 


would be 


1.0724 + @©.0391 + 9.1439 = 1.2554 kilowatts 


and the per cent of heat loss would be 


&) 
iN 


net power/overall cycle efficiency 
= 2.079/0.3389 = 6.1345 kilowatts (heater) 
% heat loss = 1.2554/6.1345 = 20.46% 


and overall cycle efficiency accounting for heat loss 


to the ocean (net efficiency) was 


= 9.3389/(1 + 0.2046) = 9.281 
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